Urban tranport phenomena in the street 

canyon 



Maciej M. Duras * 

Institute of Physics, Cracow University of Technology, ulica Podchorqzych 1, PL-30-084 Cracow, 

Poland 

"The Proceedings of the Comphy02 Computational Physics of Transport and Interface Dynamics 
conference", M. Schreckenberg et al Eds., Springer-Verlag, Berlin, Germany (2002). 

(AD 2002, 13th March) 



Abstract 

A field deterministic model of the vehicular dynamics in a generic urban 
street canyon with two neighboring canyons is considered. The assumed hy- 
drodynamical model of vehicular movement is coupled to the gasdynamical 
model of the air and emitted pollutants. The vehicles are assumed to move 
on distinct left lanes and right ones. At the upstream and downstream ends 
there are coordinated traffic lights which introduce control parameters to the 
field model. The model of optimal control of street canyon dynamics is based 
on the two optimal multi-criteria control problems. The problems consist of 
minimization of the dimensionless functionals of the cumulative total travel 
time, global emissions of pollutants, and global concentrations of pollutants, 
both in the studied street canyon, as well as together with the two nearest 
neighbor substitute canyons, respectively. 
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I. THE FIELD MODELS OF VEHICLES AND POLLUTANTS. 



In the present article we consider two coupled deterministic field model of the urban street 
canyon. The vehicular fields are one-dimensional in spatial variable x and they depend on 
time variable t. The considered fields are the vehicular number density ki vt (x,t), vehicu- 
lar velocity wf vt (x, t), emissivity of exhaust gases ef ctvt (x,t), and vehicular heat emissivity 
af vt (x,t), where vt is emission type's number, ct is number of the constituent of emitted 
exhaust gases. It is assumed that there are rii = n L left lanes s = 1 and / is the left 
lane's number, / = 1, ...,n L , whereas for n 2 = ur right lanes s = 2 and I is the right lane's 
number, I = l,...,n R , vt = 1,...,VT, ct = 1,...,CT. VT is total number of types of ve- 
hicular emissions, and CT is total number of emitted exhaust gases. The abovementioned 
model consists of the equations of balances of vehicular numbers [5] and vehicular equa- 
tions of state (Greenshields' equilibrium speed-density u-k model [6]) as well as exhaust gas 
emissivities and heat emissivities. The vehicular flow in the canyon is multilane bidirec- 
tional one-level rectilinear, and it is considered with two systems of coordinated signalized 
junctions at the upstream and downstream ends of lanes [1], [2], [3], [4]. The considered 
vehicles belong to distinct vehicular classes: passenger cars, and trucks. Emissions from the 
vehicles are based on technical measurements of the following types of pollutants: carbon 
monoxide CO, hydrocarbons HC, nitrogen oxides NO x . The field model of pollutants is 
gasdynamical. The gases are Newtonian viscous perfect and noninteracting. The consid- 
ered fields are three-dimensional in spatial variables (x, y, z) and one-dimensional in tempo- 
ral variable t and they are mixture's density p(x,y,z,t), velocity v(x,y,z,t), temperature 
T(x,y,z,t), pressure p(x,y,z,t), as well as constituents' mass concentrations Ci(x,y,z,t), 
and partial pressures Pi(x,y, z,t),i = 1,...,N, where N is total number of mixture's ele- 
ments. (N = N E — 1 + N A ). The first N E — 1 = 3 gases are the exhaust gases emitted 
by vehicle engines during combustion CO, CH, NO x . The remaining Na = 9 gases are the 
constituents of air: O2, N2, Ar, CO2, Ne, He, Kr, Xe, H2. The field equations are the balances 
of mixture's mass, linear momentum, energy, and equation of state (Clapeyron's law) as 
well as the balances of masses of constituents and constituents' equations of states (Dal- 
ton's law). The studied sources are mixture mass source S(x,y, z,t), mixture energy source 
a(x,y, z,t), and constituents' mass sources Sf(x,y,z,t) (exhaust gases and consumed oxy- 
gen). The boundary- initial problems are of mixed types (Dirichlet's and von Neumann's 
types). The equations of dynamics are solved by the finite difference scheme. The two sepa- 
rate multi-criteria optimization problems are posed by defining the dimensionless functionals 
of cumulative total travel time, global emissions of pollutants, and global concentrations of 
pollutants, either in the studied street canyon or in the canyon and its two nearest neighbor 
substitute canyons. The vector of control is a five-tuple composed of two cycle times, two 
green times, and one offset time between the traffic lights. The optimal control problems 
consists of minimization of the two functionals over the admissible control domain manifold. 



II. THE GOVERNING EQUATIONS. 



Under the abovementioned specifications, the set of governing equations is formulated 
as follows [7], [8]: 
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where wf^j is maximum free flow speed, kf vt - iaia is jam vehicular density, rj is the first 
viscosity coefficient for the air, £ is the second viscosity coefficient, F = pg is the gravitational 
body force density, g is the gravitational acceleration, D im = D mi is the mutual diffusivity 
coefficient from the z-th constituent to m-th one, and Du is the autodiffusivity coefficient of 
the i-th constituent, and kT, m is the thermodiffusion ratio of the m-th constituent, e is the 
mass density of intrinsic (internal) energy of the air mixture, T is the stress tensor, symbol 
: denotes the contraction operation, q is the vector of flux of heat. We assume that [2]: 
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where €j is the mass density of intrinsic (internal) energy of the ith constituent of the 
air mixture, rrii the molecular mass of the zth constituent, ks is Boltzmann's constant, 
Hi is the complete partial chemical potential of the ith constituent of the air mixture (it 
is complete since it is composed of chemical potential without external force field and of 
external potential), m air = 28.966 [u] is the molecular mass of air, S m k is Kronecker's delta, 
c Pt i is the specific heat at constant pressure of the ith constituent of air mixture, h is Planck's 
constant, jj is the vector of flux of mass of the ith constituent of the air mixture, and k is 
the coefficient of thermal conductivity of air. These magnitudes were derived from Grand 
Canonical ensemble with external gravitational Newtonian field. 



III. OPTIMAL CONTROL PROBLEMS. 

Let us consider the measures of the total travel time (TTT): J TTT [5], emissions (E): J E , 
and concentrations (C): Jc of exhaust gases in the street canyon as well as in the canyon and 
two its nearest neighbor substitute canyons, </TTT,ext </E,ext <^c,ext, respectively. Therefore the 
appropriate optimization problems may be formulated as follows [2]. The measures depend 
on traffic lights at the upstream and downstream ends of the lanes. The vector of boundary 
control u reads: 

u=(g 1 ,C 1 ,g 2 ,C 2 ,F)EU adm , (18) 

g m are green times, C m are cycle times, F is offset time, and [f adm is a set of admissible 
control variables. The TTT, E, C functionals are given by: 
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The integrands fcf t , e. s l ct vt , c* in the functional depend on the control vector u through the 
boundary conditions, through the equations of dynamics, as well as, through the sources. 
The value of the vector of control u directly affects the boundary conditions for vehicular 
densities, velocities, and emissivities. It also affects the sources. Next, it propagates to the 
equations of dynamics and then it influences the values of functionals. Pstp is the density of 
air at standard temperature and pressure STP, c^stp is concentration of the zth constituent 
of air at standard temperature and pressure, a, b, c, are dimensions of the street canyon, T$ 
is time of simulation. 

We define two additional functionals. The first one is the composite total travel time, 
emissions and concentrations of the pollutants in the single canyon, whereas the second one 
in the canyon and two its neighbors: 

H(sl, u) = Q!ttt-Pttt(u) + a E F E (u) + a c F c (u). (25) 
tf cxt (a,u) = tf(a,u)+ (26) 

+aTTT,ext«/TTT,ext(u) + OlE.ext ^E.ext (u) + ac,cxt«/c,ext(u). 

The scaling parameters «ttt, «e, &c, «TTT,ext, «E,ext, «c,ext make the two functionals dimen- 
sionless, and 

a = (a TT T, a E , «c) e ^ adm , (27) 
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where A adm , A a xt m are the sets of admissible scaling parameters. We formulate two separate 
multi-criteria optimization problems consisting in minimization of the functionals H, H cxt 
with respect to control vector u over admissible domain and with respect to scaling param- 
eters, while the equations of dynamics are fulfilled: 

H* = H(sl\ u*) = min{a G A adm , u G £/ adm : if (a, u)}. (29) 



ff e * xt = // cx t(a: xt , < xt ) = min{a cxt G u G £/ adm : ff ext (a, u)}, (30) 

where H*, ff* xt , are the minimal values of the two functionals, (a*, u*), and (a* xt , u* xt ), are 
scaling and control vectors at which the functionals reach the minima, respectively. 
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V. CONCLUSIONS. 

The proecological urban traffic control idea and the model of the street canyon have been 
developed. It is assumed that the proposed model represents the main features of complex 
air pollution phenomena. 
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